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We present a theory of the metal-insulator transition in a disordered two- 
dimensional electron gas. A quantum critical point, separating the metallic 
phase which is stabilized by electronic interactions, from the insulating phase 
where disorder prevails over the electronic interactions, has been identified. 
The existence of the quantum critical point leads to a divergence in the den- 
sity of states of the underlying collective modes at the transition, causing the 
thermodynamic properties to behave critically as the transition is approached. 
We show that the interplay of electron-electron interactions and disorder can 
explain the observed transport properties and the anomalous enhancement of 
the spin susceptibility near the metal-insulator transition. 
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Measurements in high-mobility two-dimensional (2D) semiconductors have made it possi- 
ble to study the properties of the 2D electron gas at low carrier densities. In a clean system at 
very low carrier densities, the electrons are expected to solidify into a Wigner crystal. In a wide 
range of low electron concentrations, prior to the Wigner crystal phase, the electron system ex- 
ists as a strongly interacting liquid (1 ). Understanding the properties of this strongly interacting 
system in the presence of disorder has proven to be an extremely challenging theoretical prob- 
lem. The recent experimental discovery of a metal-insulator transition (MIT) in 2D (2), when 
none was expected for more than a decade ©, generated renewed interest in this field (4). 

The MIT, observed initially in high-mobility silicon metal-oxide-semiconductor field-effect 
transistors (Si-MOSFETs) and subsequently in a host of other 2D systems like GaAs het- 
erostructures, occurs when the resistance R is of the order of the quantum resistance h/e 2 , 
emphasizing the importance of quantum effects. In the metallic phase the resistance drops no- 
ticeably as the temperature is lowered. This drop is suppressed when an in-plane magnetic field 
is applied (4). Additionally, Si-MOSFET samples show a strong enhancement of the spin sus- 
ceptibility as the MIT is approached There are indications that the spin susceptibility in 
samples of different mobilities behaves critically near the transition. The sensitivity to in-plane 
magnetic fields together with the anomalous behavior of the spin susceptibility highlight the 
importance of electron-electron (e-e) interactions in this phenomenon. 

In the presence of impurities, perturbations of the charge and spin densities (if spin is con- 
served) relax diffusively at low frequencies and large distances. Formally, this relaxation occurs 
via the propagation of particle-hole pair modes. In a system which obeys time-reversal sym- 
metry, the modes in the particle-particle channel (i.e., the Cooper channel) also have a diffusive 
form, leading to the celebrated weak-localization corrections to the conductivity. These two 
mode families are, in conventional terminology, referred to as Diffusons and Cooperons, re- 
spectively. Taken together, they describe the low energy dynamics of a disordered electron liq- 
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uid 0. The modes have the diffusive form V(q, to) = 1/ (Dq 2 — izu>), where D is the diffusion 
coefficient, and the parameter z determines the relative scaling of the frequency (i.e., energy) 
with respect to the length scale (00); z — 1 for free electrons. The addition of e-e interactions 
results in the scattering of these diffusion modes. Both D and z therefore, acquire corrections 
in the presence of the interactions. Conversely, diffusing electrons dwell long in each other's 
vicinity, becoming more correlated at low enough energies. As a result, the e-e scattering ampli- 
tudes 7 2 and 7 C characterizing the scattering of the Diffuson and Cooperon modes, respectively, 
acquire corrections as a function of the disorder. All these corrections, due to the diffusive form 
of the propagator T>(q, ui), are logarithmically divergent in temperature in 2D (T771. They signal 
the breakdown of perturbation theory and the need for a re-summation. 

A consistent handling of these mutually coupled corrections to D and z on the one hand, 
and 72 and 7 C on the other, requires the use of renormalization group (RG) that effectively sums 
the logarithmic series 1771 172b . allowing one to approach the strong-coupling region of the 
MIT (IPll7?b. We show that an internally consistent solution of the MIT can be obtained within 
a suitably defined large- N model involving iV flavors of electrons. Bearing in mind that the 
conduction band in semiconductors often has almost degenerate regions called valleys (|7?|> , the 
electrons carry both spin and valley indices. Taking the number of valleys n v — > oo, we derive 
the relevant RG equations to two-loop order, and show that a quantum critical point (QCP) that 
describes the MIT exists in the 2D interacting disordered electron liquid. 

As the inter-valley scattering requires a large change of the momentum, we assume that the 
interactions couple electrons in different valleys but do not mix them, i.e., inter-valley scattering 
processes, including those due to the disorder, are neglected. This assumption is appropriate for 
samples with high-mobility. In this limit, the RG equations describing the evolution of the 
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resistance and the scattering amplitude 72 in 2D are known to have the form (1751) : 
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where £ = — \n(Tr/h) with Tr <C % (diffusive regime), r is the elastic scattering time and 
the dimensionless resistance parameter p = (e 2 /nh)R is related to D as p = 1/ ' (2^) 2 n v vD, 
with v being the density of states of a single spin and valley species. (For repulsive interactions, 
the scattering amplitude 7 C scales to zero when n v is finite and is therefore, neglected for now; 
the situation in the large-r^ limit is different and is discussed later.) Though the initial values 
of p and 2, determined at some initial temperature, depend on the system and are therefore not 
universal, the flow of p can be described for each n v by a universal function R(r]) {8): 



where R(r]) is a non-monotonic function with a maximum at a temperature T max corresponding 
to the point where dp/d£ in Eq. ([T]) changes sign. It follows that the full temperature dependence 
of the resistance p is completely controlled by its value p max at the maximum; there are no other 
free (or fitting) parameters. Note that the anti-localization effect of the e-e interactions funda- 
mentally alters the commonly accepted point of view (3) that in 2D all states are "eventually" 
(i.e., at T = 0) localized. 

This solution has obvious limitations, however. Because the RG equations were derived 
in the lowest order in p, the single curve solution R(r/) in Eq. © cannot be applied in the 
critical region of the MIT where p ~ 1. (In fact, for n v = 2, R(r]) describes quantitatively 
the temperature dependence of the resistance of high-mobility Si-MOSFETs in the region of p 
up to p ~ 0.5 (1751) . which is not so far from the critical region.) Therefore, to approach the 
MIT the disorder has to be treated beyond the lowest order in p, while adequately retaining 
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the effects of the interaction. This also touches upon the delicate issue of the internal consis- 
tency and nature of the theory as T — > 0. The problematic feature of the scaling Eqs. CQ) and 
© is that the amplitude 72 diverges at a finite temperature T* and thereafter the RG theory 
becomes uncontrolled (1771 17?b . Fortunately, the scale T* decreases very rapidly with n v as 
lnln(l/rT*) ~ (2n v ) 2 , making the problem of the divergence of 72 for all practical purposes 
irrelevant even for n v = 2. Still, to get an internally consistent solution up to T = it is useful 
to study the limit n v — > 00, for which T* —> 0. 

The valley degrees of freedom are akin to flavors in standard field-theoretic models. Gen- 
erally, closed loops play a special role in the diagrammatic RG analysis in the limit when the 
number of flavors iV is taken to be very large (IToT) . This is because each closed loop involves 
a sum over all the flavors, generating a large factor N per loop. It is then typical to send a 
coupling constant A to zero in the limit iV — > 00 keeping AiV finite. For interacting spin- 1/2 
electrons in the presence of n v valleys (N = 2n v ), the screening being enhanced makes the 
bare values of the electronic interaction amplitudes 72 and 7 C to scale as l/2n v . Furthermore, 
the increase in the number of conducting channels results in the resistance p to scale as l/n v . 
It is therefore natural to introduce the amplitudes 62 = 2n„72 and 9 C = 2n v ^ c together with 
the resistance parameter t = n v p, which remain finite in the large-n^ limit. The parameter t is 
the resistance per valley, t = 1/ (2n) 2 uD, and it reveals itself in the theory via the momentum 
integration involving the diffusion propagators T>(q, uj). In terms of these variables, the contri- 
bution of a closed loop connected to the rest of the diagram by one interaction amplitude, 72 (or 
7 C ), after integrating over the momentum flowing through the loop, and summing over the spin 
and valley degrees of freedom, is proportional to 2n v t^2 = t@i (or t0 c ). While, this contribution 
remains finite at large-n^, those diagrams with more than one interaction for every closed loop 
are negligible. This one-to-one correspondence between the number of loops in a given diagram 
and the number of interaction amplitudes limits the maximum number of interaction vertices 
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for a given power of t. This is the crucial simplification on taking the large-n^ limit. 

For repulsive interactions at finite-n^, re-scattering in the Cooper channel leads to the van- 
ishing of the effective amplitude 7 C at low energies. The amplitude 9 C is, however, relevant in 
the large-n„ limit as the re-scattering is not accompanied by factors of n v , We introduce the 
parameter a — 1 to mark the contributions arising from the Cooper channel. Violation of time 
reversal symmetry (e.g., by a magnetic field) suppresses the Cooperon modes ( 7). Setting a = 
switches the Cooper channel off. Following the large-r^ approximation scheme detailed above, 
the RG equations to order t 2 are derived for a = and 1: 



The constants c t = (5 - ir 2 /3)/2 « 0.8 and c e = (1 - it 2 /12)/2 » 0.08. It turns out that in 
the large-n„ limit the amplitudes 9 2 and 9 C appear together in the combination Q = 9 2 + a9 c . 
The fact that they come together as a single parameter 9 is unique to the large-r^ limit. Eq. © 
reproduces the known result to order t 2 when the electronic interactions are absent dT7b. Notice 
that while the maximum power of 6 is limited by the order of t, the opposite is not true; it is the 
number of momentum integrations of the diffusion propagators that determine the power of t. 

Eqs. © and © describe the competition between the electronic interactions and disorder in 
2D. The flow of and 0(0 are plotted in Fig.[T]for the case a = 1 (the flows are qualitatively 
the same for a = 0). The arrows indicate the direction of the flow as the temperature is lowered. 
The QCP, corresponding to the fixed point of Eqs. © and ©, is marked by the circle. The at- 
tractive ("horizontal") separatrices separate the metallic phase where t — > from the insulating 
phase where t — > oo. Crossing one of these separatrices (TT^t by changing the initial values of 
t and 9 (e.g., by changing the carrier density) leads to the MIT. It is interesting that near the 
fixed point these separatrices are almost insensitive to temperature, which is in agreement with 
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the experiments in Si-MOSFETs (Ii9ll26>l) . The accidental (but fortunate) smallness of the fixed 
point parameters t c ~ 0.3 and t c Q c ~ 0.27 permits us to believe that the two-loop equations 
derived in the large-n„ limit capture accurately enough the main features of the transition. 

We now discuss the renormalization of the parameter z due to the e-e interactions. The RG 
equation for z is described by an independent equation, which quite generally (11) takes the 
form dlnz/dt; = (3 z (t, 72, 7 C ; n v ). Observe that f3 z , as well as Eqs. © and © for t and 6, 
are all independent of z. Consequently, z in the vicinity of the MIT is critical, i.e., z ~ l/T^ 
with a critical exponent ( equal to the value of the function (3 Z at the critical point (|9llZ21)- The 
parameter z, being related to the frequency renormalization can be interpreted as the density of 
states of the diffusion modes, and therefore controls the contribution of the diffusion modes to 
the specific heat, i.e., C v ~ {zv)T Hence, in the critical regime C v ~ T 1- ^. (The fact 

that Cy must vanish as T — > constraints the exponent £ < 1.) Furthermore, the Pauli spin 
susceptibility x m the disordered electron liquid has the general form \ ~ { zu )0- +72) (QUED, 
i.e., apart from the Stoner-like enhancement of the g-factor, (1 + 72), it is proportional to the 
renormalized density of states z. The enhancement of z near the QCP implies that the spin 
susceptibility diverges as % ~ while the spin diffusion coefficient D s = D/z(l + 72) 

scales to zero. Note that at the fixed point the ^-factor remains finite and therefore the divergence 
of the spin susceptibility is a priori not related to any magnetic instability. In the large-n^ limit, 
we find that no terms of the order t 2 are generated in the equation for z: 

d\nz _ 

—— = te . (6) 

Therefore, we get for the critical exponent ( » 0.27, which is noticeably smaller than one. 

We propose that the existence of the QCP explains the anomalous enhancement of the spin 
susceptibility that has been observed near the MIT in Si-MOSFETs with different critical den- 
sities (|5]|2|). Additionally, an important consequence of this theory is that the compressibility 
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dn/d/j is regular across the transition, in agreement with the experiments in p-GaAs (I23II24I) . 

In the insulating phase the parameter t diverges at a finite scale £ c , indicating the onset of 
strong localization. Since the t 2 term in Eq. © is negative, it forces the interaction amplitude 
6 to vanish in the insulating phase. It can be shown that the parameter z is finite at the scale £ c . 
The vanishing of the interaction amplitude and a finite value of z makes the insulating phase 
to be similar to the Anderson insulator. 

On the metallic side, a state with decreasing resistance t — > and enhanced amplitude 
— > oo is stabilized by the electronic interactions as T — > in the large-n„ limit. It can be 
shown within this solution that z ~ 1/t, implying that the quasi-particle diffusion coefficient 
defined as D qp = D/z (|22|) behaves regularly in the metallic phase. Deep in the metallic phase, 
the enhancement in \ Wli l be observed only at exponentially small temperatures. This holds 
even in the case of finite-n^, because the scale T* at which the RG Eqs. ® and © become 
uncontrolled is immeasurably small even for n v = 2. Note that the question of the existence 
of the QCP at finite-n^ is distinct from the problem of the divergence of the parameter 72 at 
T = T*, as the two phenomena occur in different parts of the phase diagram. 

A description of the MIT is obtained within the two-loop approximation in the limit of a 
large number of degenerate valleys. Although the properties of the thermodynamic quantities 
in the critical region are obtained in the large-n„ limit, it captures the generic features of the 
MIT for any n v . Our solution gives a physical picture that is in qualitative agreement with the 
experimental situation. In particular, it is shown that the point of the MIT is accompanied by a 
divergence in the spin susceptibility whose origin is not related to any magnetic instability. 
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Figure 1 : The disorder-interaction, £-©, flow diagram of the 2D electron gas obtained by solving 
Eqs. © and © with the Cooper channel included (a = 1). Arrows indicate the direction of 
the flow as the temperature is lowered. The circle denotes the QCP of the MIT, and the dashed 
lines show the separatrices. 
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